Local lattice potentials and steady-state vacancies in ionic crystals 
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Basing on peculiarities of local potentials, the two principal trends in vacancy formation are 
revealed. The proximity to the threshold of local ionic stability due to a giant potential contribution 
of electronic derealization accounts for thermal anion vacancies typical of intrinsic semiconductors 
of AB type. On the other hand, the tendency towards equalizing the potential field results in a high 
concentration of structural cation vacancies observed in NiaSb and relative compounds. 
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Bulk electrostatic potential field plays a fundamen- 
tal role in crystals |Q, for it specifies both the electro- 
negativity [0, ||, H and charge transfer H, |5[ ||, f?J and 
eventually contributes to the cohesive energy there j^, 



|HJ . Described by lattice series 

IB B H3 IHl' bulk Poten- 
tials can nevertheless be defined as unique quantities sub- 
jected to invariant periodic boundary conditions jl3|, [TJJ. 
The particular event of point charge lattices is of special 
interest §§§@ || [§, [H], || [§, g || || be- 
cause it is relevant to spherically symmetric ions of which 
external potential effect amounts to that of point charges 
p7||. However, the symmetry of ionic potentials evident 

m 



cannot explain 



30 



33 



in this case [yj, [28|, 
the known fact that the steady-state concentrations of va- 
cancies of various ionic species differ even in the simplest 
case of diatomic crystals p5| . The nature of abnormal 
concentrations of cation vacancies in NiaSb and relative 
compounds |3^, |3?], |3^] is one more intriguing problem. 

In the present paper we discuss in detail that the effect 
of prevailing anion vacancies inherent in GaAs ]39l 4C , 
41, 42] as well as in structurally relative nitrides |4^, 44, 
45. 1G| can be understood if the potential contribution of 
electronic delocalization jTj| is taken into account. 

But it turns out that another treatment based on lattice 
potentials is appropriate to explain cation vacancies in 
question. 

To discuss local properties in ionic crystals, we consider 
a crystal composed of j ionic species per unit cell, with 
basis vectors hj and with total charges qj described by 
charge distributions pj (r) spherically symmetric for sim- 
plicity. According to 1131 14, M, the bulk electrostatic 



potential C/b(r) at a point r can be written as follows: 



where 



U h (r) = U C d(r) +<£to P , 

Pj (r')dr> 



top 
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the sum over i runs over the Bravais lattice specified by 
Ri, the prime on the summation sign implies missing the 



singularities of the summand, Vj is the volume of integra- 
tion appropriate to pj(r), v is the volume of the unit cell. 
It is important that the absolute convergence of interest 
can always be achieved for (|^) by introducing fictitious 
point charges, which in turn can be included into the set 
over j |L3| Note that if the distributions Pj(r) are 
non-overlapping, then (j^) describes the potential in the 
point-charge lattice, but (p|) can be represented as 



deloe 
top 1 



(4) 
(5) 



where the first term on the right-hand side of (Q) modifies 
the result of (||) so as to retrieve the potential symmetry 
in a point-charge lattice, but ^top° c determines the po- 
tential contribution of electronic delocalization in ions, 
providing that this potential shift is intrinsically negative 
and additive with respect to all ions in question [jfj], [3~4|| . 
Based on (|]), the bulk Coulomb energy takes the form 



(6) 



where the contribution of fictitious charges mentioned 
above vanishes by definition, so that j runs over all the 
actual ions in the unit cell [[13T p4| , 



1j Jv, 



(7) 



In the foregoing case of point-like charges, Uj are re- 
duced to the corresponding point-charge-lattice ionic po- 
tentials LT° jl| ||, but shifted by ^ D e p oc . 

It is significant that every ion should be bound ener- 
getically in the lattice. This claim implies that 



qj u: 



cff < o, 



(8) 



despite the fact that according to (||) and the overall neu- 
trality of the unit cell, the entire value of £b is indifferent 
to any uniform shift of the potential values. Hence, neg- 
ative ions are to stay in positive potential fields, so that 



the admissible negative value of ^top° c appears to be re- 
stricted for every actual ionic structure and the lattice 
instability should be expected otherwise. 

Note that ^top° c does not change relative values on the 
bulk potential map, but describes the binding energies 
relative to free ions, with including their own stability 
p8 . 49, 5Cj . Therefore, it can be essential in the problem 
of steady-state vacancies forming at high temperatures. 
To estimate this effect, let a vacancy be generated by the 
simple exclusion of an ion from its regular position with- 
out any relaxation of the surroundings (l^, [55], [5^, |j| . 
The total concentration of vacancies is assumed to be 
small enough, so that their influence on the bulk poten- 
tials H |I| |55| can be ignored at least while we are inter- 
ested in thermal vacancies. The energy loss appropriate 
to a vacancy of the jth species is then as follows: 



-qjUf. 



(9) 
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FIG. 1: The ratio r\ s / Ra b of a core ionic radius to a crystal 
one for the As 3- ion versus that for the Ga 3+ ion so as to 
result in the value of <£top° c = —3.757 n/a expected for GaAs. 
Curves 1, 2, 3, and 4 correspond to Rgb,/Ras = 0.225, 0.25, 
0.279, and 0.3, respectively. The thin solid line exhibits the 
case of r As /J?As = r Ga /-RGa- 



Based on (J9J) , the steady-state concentrations rij of va- 
cant states can be defined by the Boltzmann relation as 



Wj exp 



^vac 

v ktr)> 



(10) 



radii r Ga and r As 
tribute 



tribute to ^ oc . 



charged uniformly |l{j, ^(j, con- 
According to (0), it is easy to get 



where Wj is the number of energetically equivalent ions 
of the jth species in the unit cell, /cb is the Boltzmann 
constant, T is the temperature. The inequality EJ ac > 
follows from (|^)-([l0|) as a necessary condition for nj <C 1. 

To make the further discussion complete, we start from 
the simplest case of AB compounds described by point- 
charge lattices |47j]. Then there are only two charge 
species with q + = — q_ = k, w± = 1, and 



El ac = k(U[ 



T<t C )> 



(11) 



where U°_ is the point-charge-lattice potential on q_ . Ac- 
cording to (10) and (O), for the ratio n_/n + we derive 



n - ( 
— = expl - 



knT 



> 1. 



(12) 
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Typical energies associated with (|Tc]) arc of about 
10 5 K, so that all the values of rij are usually negligible. 



the contribution of ^top°° can shift a minimal 
ff to threshold 



However 

positive potential to threshold (|8|) and thus the value 
of n_ becomes appreciable. Presumably, this is the case 
in semiconductors such as GaAs. To evaluate this partic- 
ular event, we adopt the lattice spacing a = 0.56537 nm 
p6| and n = 3e, where e is the charge of proton. On uti- 
lizing the potential values for the ZnS structure from p4[ , 
the expected magnitude of n_ w 10~ 3 at T w 1000 K ]4l| ] 
appears if we invoke ^^p 00 = — 3.757 n/a. 

To estimate this value in terms of structural ionic pa- 
rameters, we may assume that there is a close-packed 
sphalerite structure composed of spherical ions of radii 
i?Ga and R~As, but their core regions (2^, |57|, |5l, |5jl of 
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FIG. 2: Gaussian charge distribution p(r) — exp(— r 2 ) (solid 
curve) and an equivalent uniform distribution (the hatched 
area) restricted by r = ri on , whereas the vertical heavy line 
corresponds to the crystal radius i?i on = ri oll /0.7. 



TABLE I: Local point-charge-lattice potentials, in units of 
n/a, describing M3SD within the DO3 structure and the cor- 
responding Bl potentials. Sb 3_ ions are on (3 sites in the 
traditional nomenclature p7[ , but Ni + ions occupy a and 7 
sites specified in units of a in the parentheses. 



U 



Structure 


U° (0,0,0) 


U° (- - -) 


H° (±,0,0) 


D0 3 (B1F3) 


11.06098140 


-4.07072302 


-2.91953536 


Bl (NaCl) 


3.49512919 





-3.49512919 


and ions in 


terms of Gaussian orbitals ||l], |(E 


|] we have 




V ^3/2 
Aon 7T ) 


7T 5 / 2 

^dcloc 


(14) 



The foregoing ionic parameters specified by (Jl3) and (|T|) 
and imposed on this charge distribution in Fig. 2 for com- 
parison are in favor of the above assessment. 

Note that here we discuss the tendency towards arising 
vacancies in a perfect system, but not their final states 
Both the classical problem [|| || Q of 
the n-type conductivity in ZnO and that of vacancies in 
nitrides jl3|, ^d|, [y| can be comprehended in the same 
fashion. The trends in complex defects associated with 
anion vacancies 

in mil 

arc clarified as well. 

Owing to (|§|) , the above mechanism enforces the stabil- 
ity of cation states and so prevents from cation vacancies 
as a thermal effect. It means that another nature brings 
about cation vacancies observed in a large concentration 
in NisSb j}?]]. Of course, the experimental existence of 
such a state implies the possibility of its thermodynamic 
description with a large number of relevant parameters 
p8| , but the underlying physical motive remains obscure 
within such an approach. To elucidate the reason pro- 
moting the appearance of this unusual state, we resort to 
the potential effect associated with structural features. 

The stoichiometric Ni3Sb compound has the DO3 struc- 
ture m 1^, m @ that is J ust the BiF 3 ° ne m> with 

the local point-charge-lattice potentials listed in Table || 
[ PH |34| , where the lattice parameter a f» 0.59 nm |37[ and 
the charge k on the Ni + ion are used as basic units. So, 
the a states of Ni + ions are formally much more stable 
than the 7 ones. Notwithstanding, we must propose that 
there are structural vacancies on the a sites, in accord 
with experiment. The change in the local potentials can 
then be written in the form 



C/9 = (l-c)[/9(D0 3 ) + cC7j> (Bl), 



(15) 



where c describes the effective admixture of the Bl-lattice 
potentials compiled in Table [j] as well. Keeping in mind 
that the Bl lattice is free from charges on a sites, there 
are two a sites per unit cell with the occupation prob- 
ability 1 — c on each. All the other sites are occupied 
completely, with the charge k(2c — 3) on a Sb 3 ~ ion for 
neutrality. The changes in U° and U° defined by 



-2.5 
-3 

-3.5 
-4 

-4.5 



7 


a 





0.1 



0.2 



0.3 



0.4 



FIG. 3: Local point-charge-lattice potentials [7° (line a) and 
(7° (line 7) measured in units of n/a versus the deviation c 
from the total occupation of a sites. The degenerate case 
pointed out by the dotted vertical line at c = 0.248 corre- 
sponds to the Ni7i.sSb28.s composition. 



with the potentials from Table | exhibit the principal 
tendency towards the equalization of the potentials op- 
erating on structurally nonequivalent Ni + ions, as shown 
in Fig. ||. It is surprising that the potential degeneracy 
obtained here takes place just at the chemical composi- 
tion which is the most prominent in experiments |37j. 

It is evident that apart from a large number of va- 
cant states, equal potentials on a and 7 sites are capable 
of promoting the extraordinary high Ni diffusivity J67| . 
Basing on this potential equality, the existence of such 
a state may presumably be associated with the creation 
of metallic bonds in the Ni subsystem that is much more 
effective if all three interpenetrating fee sublattices com- 
posed of Ni + ions contribute evenly to their formation. 
A large number of metallic bonds can in turn depress 
the effect of local vacancies on the conduction electrons. 
As a result, the charge transfer determining the actual 
value of k is expected to be the same for all Ni + ions, in 
support of the uniform description developed above. 

As anticipated, the same trend is inherent in metallic 
alloys of relative structure |3(| [jf| pgj] . The fact that 
the Coulomb contribution is to be less pronounced there 
may be regarded as one more evidence in favor of the 
potential equalization discussed. 

It is worth noting that the present consideration of the 
effect of relative potentials is nonetheless based on ab- 
solute values of local potentials in point-charge lattices 
pm so that their combination in form (|l5|) is infor- 
mative. On the other hand, keeping in mind that the 
value of UpCDOs) is, in absolute units, close to f7^ 8 , but 
the number of cations is larger, the effect of ^top°° de- 
fined like ( |l3|) violates relation (||) for anion states and 
so renders the stoichiometric NisSb compound unstable, 
in agreement with experiments [p7| . Dynamical effects 
associated with Ni diffusivity and contributing to ^top° c 
so as to restore (||) at finite c will be discussed elsewhere. 

I am grateful to Professor V. L. Ginzburg, Professor 
E. G. Maksimov and Professor V. G. Vaks for their en- 



4 



couragement. 



[1] 

[2] 

[3] 

[4] 
[5] 
[6] 

[7] 
[8] 

[9] 

[10] 
[11] 

[12] 

[13] 

[14] 
[15] 
[16] 
[17] 
[18] 

[19] 
[20] 
[21] 
[22] 
[23] 

[24] 

[25] 

[26] 
[27] 

[28] 
[29] 
[30] 
[31] 
[32] 

[33] 



Electronic address: |tholopov@ casper . che . nsk . sr 
M. A. Spackman and R. F. Stewart, in Chemical Appli- 
cations of Atomic and Molecular Electrostatic Potentials, 
edited by P. Politzer and D. G. Truhlar (Plenum Press, 
New York, 1981), pp. 407-425. 

R. J. Boyd and G. E. Markus, J. Chem. Phys. 75, 5385 
(1981). 

K. A. Van Genechten, W. J. Mortier, and P. Geerlings, 

J. Chem. Phys. 86, 5063 (1987). 

J. C. Wang, J. Chem. Phys. 73, 5786 (1980). 

M. Satoh and T. Taki, Phys. Rev. B 23, 6732 (1981). 

R. J. Cole and P. Weightman, J. Phys.: Condens. Matter 

9, 5609 (1997). 

A. F. Kohan and G. Ceder, Phys. Rev. B 57, 3838 (1998). 
T. C. Waddington, in Advances in Inorganic Chem- 
istry and Radiochemistry, edited by H. J. Emeleus and 
A. G. Sharpe (Academic Press, New York, 1959), Vol. 1, 
pp. 157-221. 

M. P. Tosi, in Solid State Physics, edited by F. Seitz and 

D. Turnbull (Academic Press, New York, 1964) Vol. 16, 
pp. 1-120. 

M. G. Trefry, E. N. Maslen, and M. A. Spackman, J. 
Phys. C 20, 19 (1987). 

F. E. Harris, in Theoretical Chemistry: Advances and 
Perspectives, edited by H. Eyring and D. Henderson 
(Academic Press, New York, 1975), Vol. 1, pp. 147-218. 
M. L. Glasser and I. J. Zucker, in Theoretical Chemistry: 
Advances and Perspectives, edited by H. Eyring and D. 
Henderson (Academic Press, New York, 1980), Vol. 5, 
pp. 67-139. 

E. V. Kholopov, Preprint 2001-01 (Inst, of Inorg. Chem., 
Novosibirsk, 2001). 



E.V. Kholopov, |cond-mat/020353C . 

E. Madelung, Phys. Z. 19, 524 (1918). 
H. M. Evjen, Phys. Rev. 39, 675 (1932). 

F. Bertaut, J. Phys. Radium 13, 499 (1952). 

F. E. Harris and H. J. Monkhorst, Phys. Rev. B 2, 4400 
(1970). 

E. F. Bertaut, J. Phys. Chem. Solids 39, 97 (1978). 

P. Herzig, Chem. Phys. Lett. 68, 207 (1979). 

D. M. Heyes, J. Chem. Phys. 74, 1924 (1981). 

V. Massidda, Physica B 114, 327 (1982). 

L. Reining and R. Del Sole, Phys. Status Solidi B 162, 

K37 (1990). 

A. Gonis, E. C. Sowa, and P. A. Sterne, Phys. Rev. Lett. 
66, 2207 (1991). 

D. N. Argyriou and C. J. Howard, Aust. J. Phys. 45, 239 
(1992). 

L. Vitos and J. Kollar, Phys. Rev. B 51, 4074 (1995). 

J. D. Jackson, Classical Electrodynamics (John Wiley 

and Sons, New York, 1962). 

S. K. Roy, Can. J. Phys. 32, 509 (1954). 

J. Birman, Phys. Rev. 97, 897 (1955). 

J. D. Levine and P. Mark, Phys. Rev. 144, 751 (1966). 

W. Van Gool and A. G. Piken, J. Mater. Sci. 4, 95 (1969). 

J. V. Calara and J. D. Miller, J. Chem. Phys. 65, 843 

(1976). 

H. Xiaoguang, H. Meichun, and L. Hong, J. Phys.: Con- 



[34 
[35 

[36 

[37 

[38; 

[39 

[40 

[41 

[42 

[43 

[44 

[45 

[46 

[47 
[48 
[49 

[50 

[51 
[52 

[53 

[54 

[55 

[56 

[57 
[58 

[59 

[60 

[61 

[62 

[63 

[64 

[65 

[66 

[67 



dens. Matter 4, 5977 f!992l. 

E. V. Kholopov, |cond-mat/020739e . 

W. Van Gool and A. G. Piken, J. Mater. Sci. 4, 105 

(1969). 

I. M. Robertson, J. Phys.: Condens. Matter 3, 8181 
(1991). 

0. G. Randl, G. Vogl, M. Kaisermayr, W. Biihrer, J. Pan- 
netier, and W. Petry, J. Phys.: Condens. Matter 8, 7689 
(1996). 

G. Bester, B. Meyer, and M. Fahnle, Phys. Rev. B 57, 
R11019 (1998). 

G. A. Baraff and M. Schliiter, Phys. Rev. B 33, 7346 
(1986). 

J. E. Northrup and S. B. Zhang, Phys. Rev. B 50, 4962 

(1994) . 

W. A. Oates, G. Eriksson, and H. Wenzl, J. Alloys Comp. 
220, 48 (1995). 

A. Janotti, A. Fazzio, P. Piquini, and R. Mota, Phys. 
Rev. B 56, 13073 (1997). 

P. Boguslawski, E. L. Briggs, and J. Bernholc, Phys. Rev. 
B 51, 17255 (1995). 

P. Piquini, R. Mota, T. M. Schmidt, and A. Fazzio, Phys. 
Rev. B 56, 3556 (1997). 

1. Gorczyca, A. Svane, and N. E. Christensen, Phys. Rev. 
B 60, 8147 (1999). 

W. Orellana and H. Chacham, Phys. Rev. B 63, 125205 
(2001). 

E. V. Kholopov, Zh. Struct. Chim. 43, 600 (2002). 
R. E. Watson, Phys. Rev. Ill, 1108 (1958). 
K. Doll, M. Dolg, P. Fulde, and H. Stoll, Phys. Rev. B 
52, 4842 (1995). 

E. V. Stefanovich and T. N. Truong, J. Phys. Chem. 102, 
3018 (1998). 

M. Leslie and M. J. Gillan, J. Phys. C 18, 973 (1985). 
S. Bhowmick, D. Roy, and R. Bhattacharya, Chem. Phys. 
Lett. 148, 317 (1988). 

A. M. Stoneham, V. T. B. Torres, P. M. Masri, and 

H. R. Schober, Philos. Mag. A 58, 93 (1988). 

J. R. Macdonald, A. P. Lehnen, and D. R. Franceschetti, 

J. Phys. Chem. Solids 43, 39 (1982). 

C. Wolverton and A. Zunger, Phys. Rev. B 51, 6876 

(1995) . 

R. W. G. Wyckoff, Crystal Structures, Vols. 1, 2 (Inter- 
science Publishers, New York, 1963, 1964). 
J. Callaway and M.L. Glasser, Phys. Rev. 112, 73 (1958). 
W. W. Simmons and C. P. Slichter, Phys. Rev. 121, 1580 
(1961). 

S. Olszewski and T. Kwiatkowski, J. Mol. Struc. (Theo- 
chem) 235, 391 (1991). 

L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, 1960). 

J. R. Chelikowsky and J. Burdett, Phys. Rev. Lett. 56, 
961 (1986). 

E. V. R. de Castro and F. E. Jorge, J. Chem. Phys. 108, 
5225 (1998). 

A. F. Kohan, G. Ceder, D. Morgan, and C. G. Van de 
Walle, Phys. Rev. B 61, 15019 (2000). 

S. B. Zhang, S.-H. Wei, and A. Zunger, Phys. Rev. B 63, 
075205 (2001). 

C. C. Ling, S. Fung, C. D. Beling, and W. Huimin, Phys. 
Rev. B 64, 075201 (2001). 

B. V.Reddy, D. H. Sastry, S. C. Deevi, and S. N. Khanna, 
Phys. Rev. B 64, 224419 (2001). 

Th. Heumann and H. Stiier, Phys. Status Solidi 15, 95 



(1966). Rev. B 59, 9239 (1999). 

L. Manosa, J. Zarestky, M. Bullock, and C. Stassis, Phys. 



